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Abstract 

We consider the group of permutations of the vertices of a lattice. A random walk is 
generated by unit steps that each interchange two nearest neighbor vertices of the lattice. 
We study the heat equation on the permutation group, using the Laplacian associated 
to the random walk. At t = we take as initial conditions a probability distribution 
concentrated at the identity. A natural conjecture for the probability distribution at 
long times is that it is 'approximately' a product of Gaussian distributions for each 
vertex. That is, each vertex diffuses independently of the others. We obtain some formal 
asymptotic results in this direction. The problem arises in certain ways of treating the 
Heisenberg model in statistical mechanics. 
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This paper represents work in progress, and is written in a casual manner, as notes 
for a seminar or perhaps a physics article. But the work is interesting and non-trivial, 
and perhaps will inspire research, many natural questions will appear. We of course will 
be clear about what is proved and what is conjectured, some proofs will be sketched. 
As will be seen this work is somewhat discouraging about the application sought to the 
Heinsenberg model problem, but opens some pleasant mathematical vistas. 

We study a d-dimensional periodic lattice cube, A of edge size L. The number of its 
vertices, V, is given as 

|A| = N = L d = ^{V}. (1) 

We set Q to be the permutation group on V. H is the element in the group algebra given 

as 

H = -J2(h 3 -I). (2) 
in- 
here % ~ j indicates that % and j are neighboring vertices in the lattice. I is the identity 
element in Q, and is the group element that interchanges vertices % and j leaving the 
other vertices alone. H is thus the Heisenberg model Hamiltonian "promoted" from an 
operator in the Hilbert space to an element of the group algebra. (One need not be 
familiar with the Heisenberg model for purposes of this paper.) 

We then consider the group algebra element e~ Ht . This can be expanded as a linear 
combination of group elements, the g p 

e" m = E f(9 P ,t)g P (3) 

v 

where p labels the group elements. Each g p represents a mapping of the vertices, 

g p :i~^ p(i) (4) 

so p is specified by 



(P(1),P(2), ) (5) 



Equation (5) is a point in (Z d ) N , in fact in (A)^. The cardinality of the set of such 
points is N\. Such points in (A)^ are restricted by the condition that all the p(i) in (5) 
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are distinct. This subset we label (A)^*. The respective cardinalities of (A)^* and (A)^ 
are AH and N N . The latter space is simpler, being a periodic lattice cube. 

We have a time dependent probability function f(g,t) on the permutation group (or 
on (A)^*). We now extend f(g,t) to f e (x,t),x in (A)^. Of course such extensions are 
not unique. The motivation for extending / will soon be clear, largely depending on the 
simplicity of (A)^ over (A) N *. The extended function will no longer be a probability 
function. 

The relation between the Heisenberg model of ferromagnetism and our random walk 
on the permutation group was beautifully developed by R. Powers in [1]. Inspired by 
this work, the author presented a possible avenue towards proving the phase transition 
of the Heisenberg model in [2], a development of ideas in [1]. The central relation needed 
in the proof envisioned in [2] is of the form 



A is the lattice Laplacian on A. (Of course, equation (6) need only hold on f(g,t), 
but we presume the extension satisfies (6).) We do not now make explicit the degree of 
approximation implied by =. The right side of (6) is a product of gaussians (associated 
to independent random walks). It also is a solution of the heat equation on (A)^ with 
its natural lattice Laplacian! 

At the very least we would want (6) to imply 



r(x,t)=C N H (e A ') w fort large. 



(6) 



lim / e (£,0=CVlimn (e At K 



(7) 



We note 



lim f(g,t) 

t— >oo 



1 



(8) 



and 



lim (e At ) 



1 



(9) 



N 



Restricting to points x in (A)^* equation (7) becomes 



From (10) and Stirling's formula we have determined the Cn in (6) satisfy 



lim (C N ) N = e 

N—>oo 



(11) 



We will find an extension f e of / satisfying a differential equation 

df e 



dt 



a r + v r 



(12) 



A is the usual lattice Laplacian on (A) N a periodic version of the lattice Laplacian on 
Z dN . Here we will only consider V built up with "two-particle" interactions. Thus 
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(13) 



where Vij describes the interaction of the two vertices that at time t = were at i and 
j respectively. The V must be such that the solution of (12) (with initial conditions, the 
identity at t — 0) restricted to (A) N * agrees with / as defined in equation (3). It indeed 
is possible to find such V that achieve this. 

A form for Vij that works is given as follows. We apply this V to a product function 

m = n <p t ( Xl ) (i4) 



{V applied to product functions determines V uniquely.) 

(v^) (x) = - n <M^)-E-£ 



x z,V Xj,y+i 1 x j,V Xi,y+i 1 XtiV x j,V 1 Xi ,y+i Xj , y+i 



My) -<f>i(y + i)] ■ \<f>j(v) - faiv + 5 



(15) 



r is arbitrary. There are d orthonormal unit vectors, i. So the paths interact only when 
the vertices are in nearest neighbor position. The expression (15) certainly is uniquely 
determined by the conditions above. This is straightforward to show.... though I labored 
weeks on it. Herein we work with r set equal to zero. (But we believe working with an 
r 7^ at least if |r| < 1, leads to no essential changes in the form of our calculations and 
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results. We have studied this r^O situation a little, one need consider more diagrams 
than in the r = case. We will briefly comment on this again later.) 

We return to equation (6) for some further deliberations. We consider summing x 
over (A) N = A and over (A) N * = B. 

f e &t) =C N J2 IKOw = C N for t large (16) 
^ f e (x,t) = £ f(x,t) = 1 = CV £ U(e A % Xi for t large (17) 

x€B xGB xeB iev 

We will want the content of (16) and (17) to be given by the following conjectures. 



Conjecture 1 

fe(£/ e (^)) =(Cn) 1/N (18) 
where is given by (10) and the limit is uniform in N. 



It is the uniformity requirement that makes the conjecture most difficult. 
Conjecture 2 



l/N 

ftsJE n(e A %j =(c N y* (19) 



yxeB iev 

where is given by (10) and the limit is uniform in N . 



We believe Conjecture 2 is not very difficult to prove, and we plan to turn to it soon. 
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We now consider the solution of (12), treating V as a perturbation in the form of a 
Rayleigh-Schro dinger expansion. 

00 ft ft ft2 

f e (x, t) = V / dt n I dt n -! ■■■ dh e^t-t^y^-t^) . . . y e A h _ ( 2Q) 

n=0 Jo Jo Jo 

The right side of (20) is naturally represented as a sum of contributions of diagrams. In 
some more detail as a sum of products of the contributions of connected diagrams. All 
the computations of this paper deal with results for the sum over final states (in (20) 
over x in (A) N ) for connected tree-graph diagrams. If in equation (6) we sum over final 
states, over x, on both sides we get 

J2f e (x,t)=C N for t large. (21) 

X 

(See Conjecture 1.) 

Yes, in this paper we study the truth of (21), a much weakened form of (6). But we 
expect that if we can get on top of (21) we are well on our way to treating (6). Studying 
the decomposition of (20) into connected diagram contributions puts us in the "cluster 
expansion" framework familiar in statistical mechanics. A very complete treatment is in 
[3], but the level of sophistication of an undergraduate course in statistical mechanics is 
more than adequate. 



Theorem 1 We consider a two-particle (vertex) connected tree-graph contribution to 
(20). The diagram is an H shaped figure. The bottom legs end at vertices i and j and 
the upper arms end at Xi and Xj. The cross segment represents an interaction at time t 1 . 
The contribution of this diagram, with Xi and Xj summed over A 2 , but before integrating 
over ti is 

rl 

(22) 



t=tl 



where 0i satisfies 

d 



dt <Pi = A0x (23) 
0i (x, 0) = 6 Xii (24) 
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with A in (23) the Laplacian in A. Similarly for 2 . 



This is a simple computation patterned on the continuum or lattice version of 



-2 V0iV0 2 = / (0iA0 2 + 2 A0!) = / (0!0 2 + 2 0x) = / -Uifo) (25) 



We will always be dealing with tree-graph diagrams. Non tree-graph diagrams have 
contributions that fall off with t (for r = 0), and we are considering the t — > oo limit. 
We do not control the uniformity of this limit for the sum over all such diagrams, one 
reason for the formal nature of our computations. Our final comment on the r / 
case is that corresponding to Theorem 1, in this case, one must consider all "ladder" 
diagrams to get the same formal estimate, and not just the single tree graph. Choos- 
ing r / leads to much more work and no gain. To do better than our "results", 
if that is possible, one must consider potentials V with other than two-body forces. 



Theorem 2 will be the analog of Theorem 1 for connected diagrams involving n par- 
ticles, n not necessarily 2. 



Theorem 2 We consider all n-particle (vertex) connected tree-graph contributions to 
(20) involving vertices z±, z 2 , z n at t = 0. We sum over final positions, over A n . We 
do not integrate over t\, and over the other times in the order t n , £ n _i, t 2 . Using 
Theorem 1, and Theorem 2 inductively on n, each of these integrals will be of an explicit 
time derivative. In evaluating these integrals we keep only the lower limit, as if in 



we keep only the —f(a) term. We will collect contributions of the upper limits later. The 
"contribution of lower limits" is 





(26) 



(-inn-i)i^E ft Mv) 



t=ti 



(27) 
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where 

d 



m & = (28) 
My, 0) = 5 y , Zi (29) 



Proof. The proof is by induction on n. Let the earliest interaction corresponding to t\ be 
between vertices n — 1 and n. (That is, one of these vertices at z n at t = and the other 
at z n -i.) Upon integrating over the later times and keeping only lower limits we have 

— * 

terms with vertices 1, 2, at y, and vertices j + 1, n — 2 at y + i, where vertices n 
and n — 1 are at these two points at t±. One sums over the value of j, the points y, the 

— * 

unit vectors i, and the permutation of different possibilities for the vertices attached to y 
and y + i. We let S stand for the sum of the (n — 1)! permutations of vertices 1, n — 1. 
The following telescopic sum relation is the heart of the proof. 



EE* 



n-2 



n-i(y) + i 




E My)-- ( f ) j(y)^j+i(y + i)--4n-2{y + 1) 

3=0 

= E E 5 [<My)-0n-i(2/) - 0i (y + i)-4n~i(y + Tj\ (My) - My + 1)) 
« r 

= E E( n - !) ! [0i(y)-0n-i(y) - My + *)-<f>n-i{y + *)] (0n(y) - 0n(y + ij) (30) 

W i 

This will represent the contribution from terms where the first interaction involves vertex 
n. One then sums over the n possibilities for this first vertex coupled. (There is a factor 
of 2 in the first, t — t±, interaction we have absorbed against the fact that we are counting 
double since either end of the interaction at t — t±, could have been called n.) Thus the 
proof is short. Even writing out the details which we have raced over. But the proof is 
tricky enough, so that it's hard to be sure you're right. Counting is hard. 



-My + i 



Theorem 1 may be included in the statement of Theorem 2 as the n = 2 case. In the 
next theorem we take the same contributions as in Theorem 2 but in addition integrate 
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over t\ from to t, and sum over Z2, z n but requiring that z±, z n be a point in (A) n *. 
We first define 

a t = (31) 



Theorem 3 Let T n (t) be all the contributions considered in Theorem 2 for given n in 
addition integrated over t\ from to t, and summed over Z2, ■■■,z n with the requirement 
that zi, z n be a point in (A) n *. 

lim T n (t) = a n _! (32) 

t— >oo 



Proof. We detail the proof for n = 2 which contains all the essential points. From 
(27) to (29) we have upon integrating (27) from to t and summing over z 2 

E (e My,t)Mv,t)-Y, MMMm) (33) 

22^zi V y y J 

with 0i (y , 0) = 5 V:Z1 , 02 (y, 0) = 8 y , Z2 . 

The second term in (33) therefore vanishes, leaving 

E J2My,t)My,t) (34) 

22^21 y 



which equals 



The first term is 1 by 



EE^'Wy-'l-E (35) 

22 2/ 2/ 

E = 1 (36) 

22 

E<M^) = L (37) 
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and the second term in (35) gets to zero with t. This concludes the proof for n — 2. 
For n — 3, say, the initial points 21,22,23 and 2i,2 3 ,2 2 lead to the same set of diagram 
contributions explaining the lost (n — 1)! factors when pursued. 



We now define quantities Ai defined recursively from the a, of (31). We set 



3=0 

a formal power series in t. Then the Ai are defined by 

A = 1 

i-l 

A i = ai + Y. a k c °ef (P 2fe+1 , f - fc ) , i > 



(38) 



(39) 
(40) 



k=i 



where using Maple notation, coef (/, t s ) picks out the coefficient of t s in the formal power 
series /. This is the procedure by which we found the Aj. Actually, with Ai — 1,A 2 — 
2, A 3 = 5, A4 = 14...., the Ai are the Catalan numbers, given as 



1 



( 2i 



(41) 



Theorem 4 LetT n (t) be the analog ofT n (t) of Theorem 3 but now including upper limits, 
the whole megillah. Again integrating over all times, summing final states over (A) n ; keep- 
ing one initial vertex fixed and summing the other vertices at t = over points lying in 
(A) n *. Then 

lim f n (t) = A n _! (42) 



On the Proof 



The computation of the right side of (42), arising as a solution of (40), was perhaps 
the most difficulty and tricky business I have ever been associated with. Also I would 



10 



find it extremely difficult to write a presentable proof. Perhaps someone can come up 
with a reasonable proof. (Skeptics may prefer to call Theorem 4 a conjecture, but it is 
certainly true.) I content myself here with some points on the computation of A 2 . 

Equation (40), for % = 2 becomes 



The contributions of all contributing diagrams,when only the lower limit of the t 2 inte- 
gration is kept, is a 2 = — 1, the first term in (43). This by Theorem 3. 

Keeping the upper limit on the t 2 integration involves us with three cases. 

1) Case 1, associated to t x is V Z1:Z2 , associated to t 2 is V Z1:Z3 . 

2) Case 2, associated to ti is V zljZ2 , associated to t 2 is V Z2:Z3 . 

3) Case 3, associated to ti is V Z2:Z3 , associated to t 2 is V Z1:Z2 . 

With this notation there is a sum over z 2 and z% with the restriction z 2 ^ z% ^ Z\ ^ z 2 . 
Here the contribution of -21,-22, -23 does not equal the contribution of zi, 23, z 2 . 

Each of these three cases contributes a factor a\ = 1 to equation (43). Case 3 is the 
most interesting, and we will deal with this one case. 

The contribution of Case 3 may be represented as 



x±, x 2 , £3 lie in A 3 and -21,-22, -23 are restricted to A 3 *. x±,x 2 , x% are positions of the vertices 
at t — t and -21,-22,-23, the positions at t — 0. Recall we are keeping only the upper limit 
in the integral over t 2 , getting 



A 2 + a 2 + 3ai. 



(43) 




(44) 




(45) 
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k is the kernel of a two- vertex diagram with a single interaction at t — t\. We rewrite 
this as 

E(e A ') 2 x fe E f dt 1 k(x 1 ,x 3 ,z 2 ,z 3 ,t 1 )) . (46) 

x\ Z1,X1 \ x 3 z 2 ,z 3 J ° / 

We wish to compare this expression to 

E( eAt ) z x { E E f *i*(*i,*3,^,*i)] • (47) 

XI Zl ' Xl \X1,X3 Z 3 J0 J 

= 1-1. (48) 

by Theorem 3. Using translation invariance of the kernel k, (46) and (47) differ by 
an error that goes to zero with t (from the different restrictions on the z's in the two 
expressions). In (47) we neglect the restriction that z 2 and z% may not equal z\. 



We turn to the relation of Conjecture 1, equation (18) 

l/N 



hm (£/ e (f,t)| ={C N f N . (49) 

\xeA 



We work in the limit N large, and t — > oo (before TV — > oo). We thus want 

(ehm)) = e ( 5 °) 



and eschew considering uniformity oft limit with respect to N. We view f e (x, t) expressed 
a sum of products of connected diagrams. For N large we expect this sum to be dominated 
by terms with some fixed number of connectivity patterns. That is, in terms kept in the 
product there are X\N two-connected terms in the product, x 2 N three-connected terms 
in the product, x 3 N four-connected terms in the product, and so on. The expression for 

^2f e (x,t) in this limit is 

xeA 
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n 



'Aii\\ XiN N ! 



N' J (Ei(i + l)iVa; l )!(iV - £(i + 1)A^ 

(E(i + l)^i)! 



(51) 



n(xiiv)! n((i + i)!) XiJV ' 

The first set of parentheses includes the contributions of the diagrams t — > oo limit, from 
Theorem 4. The next factor, a ratio of factorials, sums over which set of initial vertices 
are included in the set that are connected to other vertices. The final ratio of factorials 
sums over the connectivities of the vertices (which vertices are connected with which 
vertices) . 

Maximizing (51) over the choice of x« one finds 

lim (E r&t)) =e" (52) 



t^oo 



with 

q = -l + J2 A^--£np (53) 

where 

1 = ]T A iP i+1 . (54) 

i=0 

One wants q = 1 but (53) and (54) do not yield q — 1. (It is not clear how to define a 
solution of (54) for p, but no reasonable definition works.) At this point we have reached 
complete frustration! 

Sometime after arriving at this impasse, we decided to consider a random walk not 
involving all the lattice vertices, but rather a fraction p of the vertices, "uniformly dis- 
tributed". This is achieved by "integrating out" a fraction (1 — p) of the vertices in the 
probability function f(x,t). This is done before the extension to f e !. It is easy to make 
the corresponding changes in all the computations of this paper, a matter of a day or 
two. 
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Replacing (52), (53) and (54) one finds 

lim(£r(a;,t)) 1/Ar = e«~ (55) 

t— >oo 

with f e here depending on pN vertices and with 

q = -1 + Y A iP i+1 — inp (56) 

~k i+l ' 



and 



i = £ a p*y . (57) 

i=0 



Where before one wanted q — 1, here we desire 

? = l + ^Ml-p). (58) 
P ' 

The miracle that happens is as follows. For p < 1/2, equation (57) is satisfied with 

p = 1 - p (59) 

and substituting this expression for p into (56) one finds (56) and (58) yield the same 
formal expansion in powers of p, valid for p < 1/2. It is interesting to note that (56), 
(57) and (58) determines both (59) and the A; L (expanding p and q in powers of p). 

If we set 

/(pp) = E 4pV (60) 
i=i 

we can solve, from (57) and (59) 

p + pf(pp) = l (61) 
p = 1 - p (62) 

to get 

/(*) = ^^ZH (63) 
and see the singularity at 2; = 1/4, or p = |. 
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So we are led to believe that Equation (6) may still be true in some suitable sense, but 
that a perturbation expansion development as undertaken in this paper is not promising. 
Presumably as p increases to value 1/2 one must consider diagrams of arbitrarily high 
connectivity 

Many interesting questions suggest themselves, of which we choose two. For p < 1/2, 
control the perturbation expansion, and obtain uniformity in N of the t — > oo limits. 
Find some way of treating p > 1/2. The first question is likely a problem about which 
to develop several Ph.D. theses. The second still requires some further ideas to gauge its 
difficulty. 
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